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Abstract 

We carry out a time dependent numerical simulation where both the hydrodynamics and the radiative 
transfer are coupled together. We consider a two-component accretion flow in which the Keplerian disk is 
immersed inside an accreting low angular momentum flow (halo) around a black hole. The injected soft: 
photons fi-om the Keplerian disk are reprocessed by the electrons in the halo. We show that in presence 
of an axisymmetric soft-photon source, the spherically symmetric Bondi flow losses its symmetry and 
becomes axisymmetric. The low angular momentum flow was observed to slow down close to the axis and 
formed a centrifugal barrier which added new features into the spectrum. Using the Monte Carlo method, 
we generated the radiated spectra as functions of the accretion rates. We find that the transitions from 
a hard state to a soft state is determined by the mass accretion rates of the disk and the halo. We 
separate out the signature of the bulk motion Comptonization and discuss its significance. We study 
how the net spectrum is contributed by photons suffering different number of scatterings and spending 
different amounts of time inside the Compton cloud. We study the directional dependence of the emitted 
spectrum as well. 

1 Introduction 

The spectral and timing properties of a black hole candidate give away the most vital clues to the under- 
standing the nature of the invisible central object. The spectrum of radiation, particularly in high energies 
give information about the thermodynamic properties of matter accreting onto a black hole. The timing 
properties give information about how these thermodynamic properties are changing with time. The ther- 
modynamic properties such as the mass density, temperature etc. and the dynamic properties such as the 
velocity components are the solutions of the governing equations. Thus, a thorough knowledge of the spectral 
and timing properties are essential (e.g., Chakrabarti, 1996). 

There are several papers in the literature which have devoted themselves to study the spectral and 
timing properties of the accretion flows around black holes. Sunyaev & Titarchuk (1980) suggested that the 
explanation of the emitted spectrum requires the presence of a Comptonizing hot electron plasma along with 
the standard disk of Shakura & Sunyaev (1973). There are several models in the literature, such as the hot 
corona on a Keplerian disk (Haardt & Maraschi, 1993), unstable inner edge of the standard disk (Kobayashi 
et al. 2003), hybrid EQPAIR model (Coppi, 1992) which uses both the thermal and non-thermal electrons 
which empirically describe the nature of the possible Compton cloud. Other models include those of Wandel 
and Liang (1991); Janiuk & Czerny, 2000; Merloni & Fabian, 2001; Zdziarski et al. 2003). In the so-called 
Two Component Advective Flow (TCAF) model of Chakrabarti & Titarchuk (1995), which is based on shock 



*himadri@bose.res.in 
tsudip@bose.res.in 
tkinsuk@bose.res.in 
§chakraba@bose.res.in 



1 



solutions in a sub-Keplerian flow (Chakrabarti, 1989), it was shown that the spectral properties are direct 
consequences of variation of accretion rates of the Keplerian (disk) and sub-Keplerian (halo) components. 
Subsequently, efforts were made to explain the timing properties. An important step in this direction is the 
theoretical work of Titarchuk and Lyubarskii (1995) and Lyubarskii (1997) who showed the influence of noise 
and turbulences on the power density spectrum. Meanwhile, almost at the same time, Molteni, Sponholz 
and Chakrabarti (1996) pointed out that the resonance effects between the cooling time scale and the infall 
time scale cause the Chakrabarti shocks (C-shocks) to oscillate and cause the most important feature of 
the power density spectrum, namely, the quasi-periodic oscillations (QPOs). Molteni, Toth and Kuznetsov 
(1999) showed that these C-shocks are actually stable even when azimuthal perturbations are given, though 
a vortex was shown to rotate anchoring the shocks, causing further enhancements in QPO power densities. 
This was further expanded by Chakrabarti, Acharyya and Molteni (2004) who relaxed the constraints on 
the equatorial symmetry and found that these shocks are prone to both vertical and radial oscillations of 
similar frequencies. Thus is it generally established that the sub-Keplerian flows are responsible for both 
the spectral and timing properties of the black hole candidates. This has been corroborated by several 
observations (Smith, Heindl & Swank, 2002; Wu et al. 2002, Soria et al. 2001, Pottschmidt et al. 2006, 
Datta & Chakrabarti, 2010). 

Given that the two component flows have been found to be useful to understand the spectral and timing 
properties, it will be important to carry out the numerical simulations of radiative flows around black holes 
which also include C-shocks. So far, however, only bremsstrahlung or pseudo-Compton cooling have been 
added into the time-dependent flows (Molteni et al. 1996; Chakrabarti et al. 2004, Proga 2007, Proga et 
al. 2008). Inclusion of the full-fledged Comptonization is prohibitively complex since the Comptonization 
efficiency depends on temperature and optical depth of the surrounding flow and this would depend on 
directions and time as well. In the present paper, we make the first attempt to incorporate the time dependent 
simulation result which includes both hydrodynamics and radiative transfer. We use the low angular halo 
along with a Keplerian disk. We find how the Comptonization affects the temperature distribution of the 
flow and how this in turn affects the dynamics of the flow as well. So far, our solutions have been steady. 
We obtain the outgoing spectrum of radiation as well. 

In the next Section, we discuss the geometry of the soft photon source and the Conipton cloud in 
our Monte Carlo simulations. The variation of the thermodynamic quantities and other vital parameters 
are obtained inside the Keplerian disk and the Compton cloud which are required for the Monte Carlo 
simulations. In §3, we describe the simulation procedure and in §4, we present the results of our simulations. 
Finally, in §5, we make concluding remarks. 

2 Geometry of the electron cloud and the soft photon source 

In Figs. [IJa-b), we present cartoon diagrams of our simulation set up for (a) spherical Compton cloud (halo) 
with zero angular momentum (specific angular momentum i.e., angular momentum per unit mass, A = 0) 
and (b) rotating Compton cloud (halo) with a specific angular momentum A = 1. In the first case (a), we 
have the electron cloud within a sphere of radius Rm = 2007-3 , the Keplerian disk resides at the equatorial 
plane. The outer edge of this disk is located at Rout = SOOr^ and it extends up to the marginally stable 
orbit Rms = 3rg. At the centre of the sphere, a black hole of mass lOAf0 is located. The spherical matter 
is injected into the sphere from the radius It intercepts the soft photons emerging out of the Keplerian 
disk and reprocesses them via Compton or inverse Compton scattering. An injected photon may undergo a 
single, multiple or no scattering at all with the hot electrons in between its emergence from the Keplerian 
disk and its escape from the halo. The photons which enter the black holes are absorbed. In the second case 
(b), due to the presence of the angular momentum of the flow, the spherical symmetry of the flow is lost. 
The other parameters of the Keplerian disk and the halo remains the same as in (a). 

2.1 Distribution of temperature and density inside the Compton cloud 

A realistic accretion disk is expected to be three-dimensional. Assuming axisymmetry, we have calculated the 
flow dynamics using a flnite difference method which uses the principle of Total Variation Diminishing (TVD) 



2 




Figure 1: Schematic diagram of the geometry of om^ Monte Carlo simulations for (a) A = and for (b) A = 1. 
Zigzag trajectories and velocity vectors are typical paths followed by the photons and the velocity vectors of 
the infalling matter inside the cloud. 



to carry out hydrodynamic simulations (see, Ryu, Chakrabarti & Molteni, 1997 and references therein; Giri 
et al. 2010). At each time step, we carry out Monte Carlo simulation to obtain the cooling/heating due 
to Comptonization. We incorporate the cooling/heating of each grid while executing the next time step 
of hydrodynamic simulation. The numerical calculation for the two-dimensional flow has been carried out 
with 900 X 900 cells in a 200rg x 200rg box. We chose the units in a way that the outer boundary (Rm) is 
chosen to be unity and the matter density is normalized to become unity. We assume the black hole to be 
non-rotating and we use the pseudo-Newtonian potential — 2(,--i) (Paczyhski & Wiita, 1980) to calculate the 
flow geometry around a black hole (Here, r is in the unit of Schwarzschild radius Vg = 2GMbhl c^)- Velocities 
and angular momenta are measured in units of c, the velocity of light and rgC respectively. In Figs. HJa-b) 
we show the snapshots of the density and temperature (in keV) profiles obtained in a steady state purely 
from our hydrodynamic simulation. The density contour levels are drawn for 0.65 — 1.01 (levels increasing 
by a factor of 1.05) and 1.01 — 66.93 (successive level ratio is 1.1). The temperature contour levels are drawn 
for 16.88 — 107.8 keV (successive level ratio is 1.05). 



2.2 Properties of the Keplerian disk 



The soft photons are produced from a Keplerian disk whose inner edge has been kept fixed at the marginally 
stable orbit Rms^ while the outer edge is located at Rout (=assumed to be at 300rg in this paper). The source 
of the soft photons have a multicolor blackbody spectrum coming from a standard (Shakura & Sunyaev, 1973, 
hereafter SS73) disk. We assume the disk to be optically thick and the opacity due to free-free absorption is 
more important than the opacity due to scattering. The emission is black body type with the local surface 
temperature (SS73): 
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The total number of photons emitted from the disk surface is obtained by integrating over all frequencies 
{v) and is given by, 
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Figure 2: Density (a) and temperature (b) contours inside the spherical halo in the absence of Compton 
cooling. Here, densities are in normalized unit and temperatures are in keV. A = is chosen. See text for 
details. 



The disk between radius r to r + 6r injects dN(r) number of soft photons. 

dN{r) = 4:TTrSrH{r)nj{r), 
where, H{r) is the half height of the disk given by: 
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In the Monte Carlo simulation, we incorporated the directional effects of photons coming out of the Keplerian 
disk with the maximum number of photons emitted in the z-direction and minimum number of photons are 
generated along the plane of the disk. Thus, in the absence of photon bending effects, the disk is invisible as 
seen edge on. The position of each emerging photon is randomized using the distribution function (Eq. 3). 
In the above equations, the mass of the black hole Mhh is measured in units of the mass of the Sun (Mq), 
the disk accretion rate Mdn is in units of 10^^ gm/s. We chose Mbh = 10 in the rest of the paper. 



3 Simulation Procedure 

In a given run, we assume a Keplerian disk rate (rhd) and a sub-Keplerian halo rate {rhh). The specific 
energy (e) of the halo provides the hydrodynamic (e.g., number density of the electrons and the velocity 
distribution) and the thermal properties of matter. Since we chose the Paczyhski-Wiita (1980) potential, 
the radial velocity is not exactly unity at r = 1, the horizon. It becomes unity just outside. In order not 
to over estimate the effects of bulk motion Comptonization (Chakrabarti & Titarchuk, 1995) which is due 
to the momentum transfer of the moving electrons to the horizon, we kept the highest velocity to be 1. We 
use the absorbing boundary condition at r = 1.5 (A = case) and r = 2.5 (A = 1 case). These simplifying 
assumptions do not affect our conclusions, especially because we are studying inviscid flow and the specific 
angular momentum is constant. Photons are generated from the Keplerian disk as mentioned before and 
may be intercepted by the sub-Keplerian halo (sphere in Fig. [1^ and cylinder in Fig. [T|d). 

To begin the Monte Carlo code, we randomly generated soft photons from the Keplerian disk. The energy 
of the soft photon at radiation temperature T(r) is calculated using the Planck's distribution formula, where 
the number density of the photons {n-y{E)) having an energy E is expressed by. 
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where, b — l/kT{r) and C(3) = l^^ — 1.202, the Riemann zeta function. 

Usmg another set of random numbers we obtained the direction of the injected photon and with yet 
another random number we obtained a target optical depth at which the scattering takes place. The 
photon was followed within the electron cloud till the optical depth (t) reached t^. The increase in optical 
depth (c?t) during its traveling of a path of length dl inside the electron cloud is given by: dr = pncrdl, where 
Pn is the electron number density. 

The total scattering cross section a is given by Klein-Nishina formula: 
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where, x is given by, 

x= j7 1-/^- , (7) 

Te = e^/mc^ is the classical electron radius and m is the mass of the electron. 

We have assumed here that a photon of energy E and momentum -f is scattered by an electron of energy 

jmc^ and momentum — 7771^, with ■f — I 1 — J and fi = fJ.v. At this point, a scattering is allowed 
to take place. The photon selects an electron and the energy exchange is computed using the Compton or 
inverse Compton scattering formula. The electrons are assumed to obey relativistic Maxwell distribution 
inside the Compton cloud. The number dN(p) of Maxwellian electrons having momentum between p to 
p + dp is expressed by, 

dN{p) cx expi-ip^c^ + m^c^f''^/kT^]dp. (8) 

We take a steady state flow profile from a hydrodynamics code to start the Monte Carlo simulation. 
When a photon interacts with an electron via Compton or inverse-Compton scattering, it loses or gains 
some energy (Ai5). At each grid point, we compute IS.E. We update the energy of the flow at this grid 
by this amount and continue the hydrodynamic code with this modified energy. This, in turn, modify the 
hydrodynamic profile. Thus the Monte Carlo code for radiative transport and numerical code are coupled 
together. In case the final state is steady, the temperature of the cloud would be reduced progressively to a 
steady value from the initial state where no cooling was assumed. If the final state is oscillatory, the solution 
would settle into a state with Comptonization. 



3.1 Details of the Coupling procedure 

Once a steady state is achieved in the non-radiative hydro-code, we compute the spectrum using the Monte 
Carlo code. This is the spectrum in the first approximation. To include cooling in the coupled code, we 
follow these steps: (a) we calculate the velocity, density and temperature profiles of the electron cloud from 
the output of the hydro-code, (b) Using the Monte Carlo code we calculate the spectrum, (c) Electrons 
are cooled (heated up) by the inverse-Compton (Compton) scattering. We calculate the amount of heat loss 
(gain) by the electrons and its new temperature and energy distributions and (d) taking the new temperature 
and energy profiles as initial condition, we run the hydro-code for a period of time. Subsequently, we repeat 
the steps (a-d). In this way, we get an opportunity to see how the spectrum is modified as the iterations 
proceed. The iterations stop when two successive steps produce virtually the same temperature profile and 
the emitted spectrum. 



3.1.1 Calculation of energy reduction using Monte Carlo code: 

For Monte Carlo simulation, we divide the Keplerian disk in different annuli of width D{r) = 0.5. Each 
annulus is characterized by its central temperature T(r). The total number of photons emitted from the disk 
surface of each annulus can be calculated using Eqn. This total number comes out to be ^ ^O"^^^""^ 
for riid = 1.0. In reality, one cannot inject this much number of photons in Monte Carlo simulation because 
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of the limitation of computation time. So we replace this large number of by a low number of bundles, say, 
^compif) ^ 10^ and calculate a weightage factor 

dN{r) 

Jw = 



Clearly, from each annulus, the number of photons in a bundle will vary. This is computed exactly and used 
to compute the change of energy due to Comptonization. When this injected photon is inverse-Comptonized 
(or, Comptonized) by an electron in a volume element of size dV , we assume that fw number of photons 
has suffered similar scattering with the electrons inside the volume element dV . If the energy loss (gain) per 
electron in this scattering is Ai?, we multiply this amount by fw and distribute this loss (gain) among all 
the electrons inside that particular volume element. This is continued for all the 10'' bundles of photons and 
the revised energy distribution is obtained. 

3.1.2 Computation of the temperature distribution after cooling 

Since the hydrogen plasma considered here is ultra-relativistic (7 = | throughout the hydrodynamic sim- 
ulation), thermal energy per particle is Sfc^T where ks is Boltzmann constant, T is the temperature of 
the particle. The electrons are cooled by the invcrsc-Comptonization of the soft photons emitted from the 
Keplerian disk. The protons are cooled because of the Coulomb coupling with the electrons. Total number 
of electrons inside any box with the centre at location {ir, iz) is given by, 

dNe{ir,iz) = 4:TTrne{ir,iz)drdz, (9) 

where, ne{ir, iz) is the electron number density at {ir, iz) location, and dr and dz represent the grid size along 
r and z directions respectively. So the total thermal energy in any box is given by 3kBT{ir,iz)dNe{ir,iz) = 
12TTrkBT(ir,iz)ne{ir,iz)drdz, where T{ir,iz) is the temperature at {ir,iz) grid. We calculate the total 
energy loss (gain) AE of electrons inside the box according to what is presented above and subtract that 
amount to get the new temperature of the electrons inside that box as 

AE 

ksTnewiir^iz) = kBToid{ir,iz) - — — . (10) 

MNe[ir,iz) 

3.2 Details of the hydrodynamic simulation code 

As mentioned above, after every spell of cooling by the Monte Carlo code for a very short time step, the 
hydro-code is run without assuming cooling. This procedure is repeated. While running the hydro-code the 
following process is followed. 

To model the initial injection of matter, we consider an axisymmetric flow of gas in the pseudo-Newtonian 
gravitational field of a black hole of mass M^h located at the centre in the cylindrical coordinates [R, 6, z] . 
We assume that at infinity, the gas pressure is negligible and the energy per unit mass vanishes. We also 
assume that the gravitational field of the black hole can be described by Paczyhski & Wiita (1980), 

(j){r) 



where, r = \/R? + z^, and the Schwarzschild radius is given by, 

Tg = 2GMbh/c^- 

We also assume a polytropic equation of state for the accreting (or, outflowing) matter, P = Kp"' , where, 
P and p are the isotropic pressure and the matter density respectively, 7 is the adiabatic index (assumed to 
be constant throughout the flow, and is related to the polytropic index n by 7 = 1 -|- 1 /n) and K is related 
to the speciflc entropy of the flow s. The details of the code is described in Ryu, Molteni & Chakrabarti 
(1997) and in Giri et al. (2010). 
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Our computational box occupies one quadrant of the R-z plane with < i? < 200 and < z < 200. The 
incoming gas enters the box through the outer boundary, located at Rin = 200. We have chosen the density 
of the incoming gas — 1 for convenience since, in the absence of self-gravity and cooling, the density is 
scaled out, rendering the simulation results valid for any accretion rate. As we are considering only energy 
flows while keeping the boundary of the numerical grid at a finite distance, we need the sound speed a (i.e., 
temperature) of the flow and the incoming velocity at the boundary points. For the spherical flow with zero 
angular momentum (Bondi flow), we have taken the boundary values from standard pseudo-Bondi solution. 
We injected the matter from both the outer boundary of R and z coordinate. In order to mimic the horizon 
of the black hole at the Schwarzschild radius, we placed an absorbing inner boundary at r = 1.5rg, inside 
which all material is completely absorbed into the black hole. For the background matter (required to avoid 
division by zero) we used a stationary gas with density pbg = 10~^ and sound speed (or temperature) the 
same as that of the incoming gas. Hence the incoming matter has a pressure 10® times larger than that of 
the background matter. All the calculations were performed with 900 x 900 cells, so each grid has a size of 
0.22 in units of the Schwarzschild radius. 

All the simulations are carried out assuming a stellar mass black hole (M = IOMq). The procedures re- 
main equally valid for massive/super-massive black holes. We carry out the simulations till several thousands 
of dynamical time-scales are passed. In reality, this corresponds to a few seconds in physical units. 

4 Results and Discussions 

In Table 1, we summarize all the cases for which the simulations have been presented in this paper. In 
Column 1, various cases are marked. Columns 2 and 3 give the angular momentum (A) and the specific 
energy (e) of the flow. The Keplerian disk rate (md) and the sub-Keplerian halo rate (m/i) are listed in 
Columns 4 and 5. The number of soft photons, injected from the Keplerian disk (-/Vmi) for various disk 
rates can be found in Column 6. Column 7 lists the number of photons (Ngc) that have suffered at least 
one scattering inside the electron cloud. The number of photons (Nunsc), escaped from the cloud without 
any scattering are listed in Column 8. Columns 9 and 10 give the percentages of injected photons that have 
entered into the black hole (iVf,^) and suffered scattering (p = -^^), respectively. The cooling time (to) of 
the system is defined as the expected time for the system to lose all its thermal energy with the particular 
flow parameters (namely, rhd and rhh)- We calculate to = E/E in each time step, where, E is the total 
energy content of the system and E is the energy gain or loss by the system in that particular time step. 
We present the energy spectral index a [I{E) ^ E^"] obtained from our simulations in the last column. 



Table 1: Parameters used for the simulations and a summary of results. 



Case 


A 


e 


rhd 


rhh 








Nbh [%] 


P[%] 


to [sec] 


a 


la 





22E-4 


1 


1 


4.3E-F40 


8.7E+39 


3.5E-t-40 


0.119 


20.030 


228.3 


1.15, 0.99 


lb 





22E-4 


2 


1 


1.5E+41 


2.9E+40 


1.2E+41 


0.120 


20.023 


63.6 


1.30, 1.0 


Ic 





22E-4 


5 


1 


7.3E+41 


1.5E+41 


5.9E+41 


0.121 


19.942 


12.4 


1.40, 0.96 


Id 





22E-4 


10 


1 


2.5E-F42 


5.0E-I-41 


2.0E-I-42 


0.121 


19.816 


4.2 


1.65, 0.90 


le 





22E-4 


1 


0.5 


4.3E+40 


4.7E+39 


3.9E+40 


0.070 


10.886 


380.0 


1.57 


If 





22E-4 


1 


2 


4.3E+40 


1.5E+40 


2.8E+40 


0.230 


34.324 


118.9 


1.1 


Ig 





22E-4 


1 


5 


4.3E+40 


2.6E+40 


1.8E+40 


0.502 


59.012 


48.0 


0.7 


Ih 





22E-4 


1 


10 


4.3E+40 


3.3E+40 


l.lE+40 


0.699 


75.523 


35.1 


0.45 


2a 


1 


3E-4 


1 


1 


6.3E-F40 


1.2E-f40 


5.1E-I-40 


0.285 


19.199 


79.7 


0.88 


2b 


1 


3E-4 


2 


1 


2.1E+41 


4.1E+40 


1.7E+41 


0.283 


19.278 


21.9 


0.94 


2c 


1 


3E-4 


5 


1 


l.OE+42 


1.9E+41 


8.1E+41 


0.283 


19.205 


4.3 


1.03 


2d 


1 


3E-4 


10 


1 


3.6E+42 


6.9E+41 


2.9E+42 


0.289 


18.941 


1.4 


1.17 


2e 


1 


3E-4 


10 


0.5 


3.6E+42 


3.9E+41 


3.2E+42 


0.190 


10.768 


1.9 


1.37 


2f 


1 


3E-4 


10 


1.5 


3.6E-F42 


9.3E-f41 


2.7E-f42 


0.372 


25.492 


1.1 


1.01 


2g 


1 


3E-4 


10 


2 


3.6E+42 


l.lE+42 


2.5E+42 


0.443 


30.758 


0.9 


0.95 


2h 


1 


3E-4 


10 


5 


3.6E+42 


1.8E+42 


1.7E+42 


0.690 


51.073 


0.7 


0.59 
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Figure 3: Changes in the density distribution in presence of coohng. A = and rrih — 1 for all the cases. 
Disk accretion rate rhd used are (a) 1, (b) 2, (c) 5 and (d) 10 respectively (Cases l(a-d) of Table 1). The 
density contours are drawn using the same contour levels as in Fig. [2^. 

4.1 Compton cloud with no angular momentum 

First we discuss the results corresponding to the Cases l(a-d) of Table 1. In Fig. [JJa-d) we present the 
changes in density distribution as the disk accretion rates are changed: — (a) 1, (b) 2, (c) 5 and (d) 
10 respectively. We notice that as the accretion rate of the disk is enhanced, the density distribution losses 
its spherical symmetry. In particular the density at a given radius is enhanced in a conical region region 
along the axis. This is due to the cooling of the matter by Compton scattering. To show this, in Fig. HJa-d) 
we show the contours of constant temperatures (marked on curves) of the same four cases. We notice that 
the temperature is reduced along the axis (where the optical depth as seen by the soft photons from the 
Keplerian disk is higher) drastically after repeated Compton scattering. 

In Fig. [5l^a-d), we show the hydrodynamic and radiative properties. In Fig. [S{a), we show the sonic 
surfaces. The lowermost curve corresponds to theoretical solution for an adiabatic flow (e.g., Chakrabarti, 
1990). Other curves from the bottom to the top are the iterative solutions for the Case Id mentioned 
above. As the disk rate is increased, the cooling increases, and consequently, the Mach number increases 
along the axis. Of course, there are other effects: The cooling causes the density to go up to remain in 
pressure equilibrium. In Fig. [SJd, the Mach number variation is shown. The lowermost curve (marked 1) is 
from the theoretical consideration. Plots 2-4 are the variation of Mach number with radial distance along 
the equatorial plane, along the diagonal and along the vertical axis respectively. In Fig. [SJ;, the average 
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Figure 4: Changes in the temperature distribution in presence of coohng. A = and rhh — 1 for ah the 
cases. Disk accretion rate is (a)l, (b)2, (c)5 and (d)10 respectively (Cases l(a-d) of Table 1). Contours 
are drawn using the same levels as in Fig. [5]d. 
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temperature of the spherical halo is plotted as a function of the iteration time until almost steady state is 
reached. The cases are marked on the curves. We note that as the injection of soft photons increases, the 
average temperature of the halo decreases drastically. In Fig. [SJl, we have plotted the energy dependence 
of the photon intensity. We find that, as we increase the disk rate, keeping the halo rate fixed, number of 
photons coming out of the cloud in a particular energy bin increases and the spectrum becomes softer. This 
is also clear from Table 1, Ninj increases with m^, increasing a. We find the signature of double slope in 
these cases. As the disk rate increases, the second slope becomes steeper. This second slope is the signature 
of bulk motion Comptonization. As increases, the cloud becomes cooler (Plot 5c) and the power-law tail 
due to the bulk motion Comptonization (Chakrabarti & Titarchuk, 1995) becomes prominent. 

In Fig. [5^ we show the variation of the energy spectrum with the increase of the halo accretion rate, 
keeping the disk rate (m^ = 1) and angular momentum of the fiow (A = 0) fixed. The injected multi-color 
blackbody spectrum supplied by the Keplerian disk is shown (solid line). The dotted, dashed, dash-dotted, 
double dot-dashed and double dash-dotted curves show the spectra for rhh — 0.5, 1, 2, 5 and 10 respectively. 
The injected multicolor blackbody spectrum supplied by the Keplerian disk is shown (solid line). The 
spectrum becomes harder for higher values of rhh = 1 as it is difficult to cool a higher density matter with 
the same number of injected soft photon. In Fig. we show the directional dependence of the spectrum, 
for A = 0, mil — 2, = 1 (Case If). The solid, dotted and dashed curves are for observing angles (a) 2°, 
(b) 45° and (c) 90° respectively. All the angles are measured with respect to the rotation axis (z-axis). As 
expected, the photons arriving along the z-axis would be dominated by the soft photons from the Keplerian 
disk while the power-law would dominate the spectrum coming edge-on. 

We now study the dependence of spectrum on the time delay between the injected photon and the 
outgoing photon. Depending on the number of scatterings suffered and the length of path traveled, different 
photons spend different times inside the Compton cloud. The energy gain or loss by any photon depends 
on this time. Fig. ^ shows the spectrum of the photons suffering different number of scatterings inside the 
cloud. Here 1, 2, 3, 4, 5 and 6 shows the spectrum for 6 different ranges of number of scatterings. Plot 1 
shows the spectrum of the photons that have escaped from the cloud without suffering any scattering. This 
spectrum is nearly the same as the injected spectrum, only difference is that it is Doppler shifted. As the 
number of scattering increases (spectrum 2, 3 and 4), the photons are more and more energized via inverse 
Compton scattering with the hot electron cloud. For scatterings more than 19, the high energy photons 
start loosing energy through Compton scattering with the relatively lower energy electrons. Components 5 
and 6 show the spectra of the photons suffering 19-28 scatterings and the photons suffering more than 28 
respectively. Here the flow parameters are: rhd — 1, trih = 10 and A = (Case Ih Table 1). 

In Fig. [SJl, we plot the spectrum emerging out of the electron cloud at four different time ranges. In the 
simulation, that the photons take 0.01 to 130 ms to come out of the system. We divide this time range into 4 
suitable bins and plot their spectrum. Case Ih of Table 1 is considered. We observe that the spectral slopes 
and intensities of the four spectra are different. As the photons spend more and more time inside the cloud, 
the spectrum gets harder (plots 1, 2 and 3). However, very high energy photons which spend maximum time 
inside the cloud lose some energy to the relatively cooler electrons before escaping from the cloud. Thus the 
spectrum 4 is actually the spectrum of Comptonized photons. 

4.2 Compton cloud with very low angular momentum 

We now turn our attention to the case where the cloud is formed by a low angular momentum fiow. In 
this case the flow is already axisymmetric and due to centrifugal force a weak shock wave, or at least a 
density wave would be formed. In Fig. [7{a-b), we show the contours of constant density (Fig. [7^) and 
temperature (Fig. (TJj) when no radiative transfer is included. Here the specific angular momentum of 
A = 1 was chosen. Density contour levels are drawn from 0.001 — 55.35 (the successive level ratio is 1.5), 
55.35 — 73.73 (successive level ratio is 1.1). Temperature contour levels are drawn from 2.3— 11.64 (successive 
level ratio is 1.5), 11.64 — 64.71 (successive level ratio is 1.1). We note that a shock has been formed which 
bends outwards away from the equatorial plane (Ryu, Chakrabarti & Molteni, 1997; Giri et al., 2010.). In 
Fig. IHta-d), we show results of inserting a Keplerian disk in the equatorial plane. The inner edge is located 
at 3 fg, the marginally stable orbit. Here, rhh = 1 and md —{&) 1, (b) 2, (c) 5 and (d) 10 respectively (Cases 



10 





Figure 5: (a) Sonic surfaces at different stages of iterations. The final curve represents the converged 
solution. The initial spherical sonic surface become prolate spheroid due to cooling by the Keplerian disk 
at the equatorial plane. Parameters are for Case Id (Table 1). (b) Mach number variation as a function of 
distance after a complete solution of the radiative flow is obtained. Plot no. 1 corresponds to the solution 
from adiabatic Bondi flow. Plots 2-4 are the solutions along the equatorial plane, the diagonal and the axis 
of the disk. Parameters are for Case Id (Table 1). (c) Variation of the average temperature of the Compton 
cloud as the iteration proceeds when the disk accretion rate is varied, jfih = 1. The solid, dotted, dashed 
and dash-dotted plots are for rhd = 1, 2, 5 and 10 respectively. Case numbers (Table 1) are marked. With 
the increase of disk rate, the temperature of the Compton cloud converges to a lower temperature, (d) 
Variation of the spectrum with the increase of disk accretion rate. Parameters are the same as in (c). With 
the increase in thd, the intensity of the spectrum increases due to the increase in Ninj (see, Table 1). The 
spectrum is softer for the higher value of rhd. Spectral slopes for each of these spectra are listed in Table 1. 
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Figure 6: (a) Variation of the spectrum with the increase of the halo accretion rate, keeping the disk rate 
{rhd = 1) and angular momentum of the flow (A — 0) fixed. The dotted, dashed, dash-dotted, double 
dot-dashed and double dash-dotted curves show the spectra for jrih = 0.5, 1, 2, 5 and 10 respectively. The 
injected multicolor blackbody spectrum supplied by the Keplerian disk is shown (solid line), (b) Directional 
dependence of the spectrum: A = 0, m/j = 2, riid = 1 are the flow parameters. The solid, dotted and dashed 
curves are for observing angles 2°, 45° and 90° respectively. All the angles are measured with respect to 
the rotation axis (z-axis). Intensity of spectra emerging from the cloud after suffering various number of 
scatterings (c) and at four different times (d) immediately after the injection of soft photons. Case Ih is 
assumed. The spectra of the photons suffering 0, 1-2, 3-6, 7-18, 19-28 and more than 29 scatterings are shown 
by the plots 1, 2, 3, 4, 5 and 6 (Fig. [6t) respectively, within the cloud. Curve Ih is the net spectrum for 
which these components are drawn. As the number of scattering increases, the photons gain more and more 
energy from the hot electron cloud through inverse Comptonization process. The spectra of the photons 
spending 0.01-1, 1-40, 40-100 and more than 100 ms time inside the electron cloud are marked by 1, 2, 3 
and 4 (Fig. |BJl) respectively. 
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Figure 7: Density (a) and temperature (b) contours inside the halo (A = 1) in the absence of Compton 
coohng. Densities are in normaUzed unit and temperatures are in keV. See, text for details. 

2(a-d) of Table 1). The densities used to draw the contours are the same as that in Fig. [7^. As the Keplerian 
disk rate is increased, the intensity of the soft photons interacting with the high optical depth (post-shock) 
region is increased. In Fig. |8jl, we observe that the conical region around the axis is considerably cooler. 
Thus, the density around the shock is enhanced. However, most importantly, with the increase in disk 
accretion rate, i.e., cooling, the shock location moves in closer to the black hole. This result has been already 
shown in the context of the bremsstrahlung cooling (Molteni et al. 1996). In Fig. [HJa-d), we present the 
corresponding temperatures. The parameters are the same as in Fig. |Sl[a-d) and the temperatures used 
to draw the contours are the same as that in Fig. [7}d. The Comptonization in the shocked region cools 
it down considerably. Otherwise, not enough visible changes in the thermodynamic variables are seen. To 
understand the detailed effects of the radiative transfer on the dynamics of the flow, we take the differences 
in the pressure and velocity at each grid point of the flow for Cases Id and 2d of Table 1. In Fig. [TUFa-b) 
we show the difference between the results of a purely hydrodynamical flow and the results by taking the 
Comptonization into account. Fig. \Wh is for the flow with no angular momentum and Fig. llOb is drawn 
for the specific angular momentum A = 1. The contours are of constant AP = Pc — Pa, where P is the 
pressure and the subscripts c and a represent the pressure with and without cooling respectively. The arrows 
represent the difference in velocity vectors in each grid. As expected, in both the cases, the changes are 
maximum near the axis. The fractional changes in pressures and velocities are anywhere between ^ (outer 
edge) and ^ 25% (inner edge and near the axis). Because of shifts of the shock location towards the axis, 
the variation of the velocity is also highest in the vicinity of the shock. Most importantly, the matter starts 
to fall back after losing the outward drive. This is the Chakrabarti-Manickam mechanism (2000) which is 
believed to decide the nature of the light curves of objects containing outflows. Thus we prove that not only 
the symmetry is lost by the insertion of an axisymmetric soft photon source, the cooling process also plays 
a major role in deciding the dynamics of the flow. 

We now turn our attention to the dynamical variables and spectral behaviour of the rotating flow. In Fig. 
Illb . we show the variation of the average temperature of the Compton cloud as a function of the iteration 
time of the coupled code (Cases 2(a-d), Table 1). With the increase in the disk rate, the temperature of 
the Compton cloud saturates at a lower temperature. Fig. [TTb shows the effect of the decrease in cloud 
temperature on the spectrum due to the increase in disk rate. As we increase the disk rate, keeping the halo 



13 




Figure 8: Change in the density contours in presence of coohng (A = 1). See, text for details. The conical 
region between the axis and shock wave becomes denser as the accretion rate of the Keplerian disk is 
increased. 
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Figure 9: Change in the temperature contours in presence of coohng. The parameters are the same as in 
Fig. |Sl[a-d). The temperature values used to draw the contours are the same as in Fig. [71d. Note that the 
shock shifts closer to the axis with the increase in disk accretion rate. 



200 
180 
160 
140 
120 
N 100 
80 
60 
40^ 
20 




'.ill t KftT? " ; s 



iiiiiii ' ' 



200 




20 40 60 80 100 120 140 160 180 200 20 40 60 80 100 120 140 160 180 200 
R R 



Figure 10: Difference in pressure and velocities between the flow with Comptonization and without Comp- 
tonization. Other parameters remain exactly the same. The cases are (a) Case Id and (b) Case 2d of Table 
1 respectively. 
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Figure 11: Variation of (a) average temperature of the Compton cloud with iteration time and (b) spectrum 
with the increase of disk accretion rate. A = 1 and rrih = 1 are used. The sohd, dotted, dashed and dash- 
dotted plots are for rhd = 1, 2, 5 and 10 respectively. With the increase in the disk rate, the temperature of 
the Compton cloud saturates at lower temperature. The solid, dotted, dashed and dash-dotted curves show 
the spectrum for = 1, 2, 5 and 10 respectively. The spectrum is softer for higher value of rhd- 

rate fixed, the spectrum becomes softer. 

In Fig. [T^ we show the effects of the increase of the electron number density (due to the increase of 
rhh) for a fixed disk rate. The spectrum becomes harder as we increase the halo rate keeping the number of 
injected soft the photons the same. 

We observe that the emerging spectrum has a bump, especially at higher accretion rates of the halo, at 
around 100 keV (e.g. the spectra marked Ig, Ih in Fig. [6^ and the spectrum marked 2h in Fig. [T2|) . A 
detailed analysis of the emerging photons having energies between 50 to 150 keV was made to see where in 
the Compton cloud these photons were produced. In Fig. [T^ . we present the number of scatterings inside 
different spherical shells within the electron cloud suffered by these photons (50 < E < 150 keV) before 
leaving the cloud. Parameters used: rhd — 1, i^h — 10 and A = 0. The light and dark shaded histograms 
are for the cloud with and without bulk velocity components, respectively. We find that the presence of 
bulk motion of the infalling electrons pushes the photons towards the hotter and denser [Figs. Ufa-b)] inner 
region of the cloud to suffer more and more scatterings. We find that the photons responsible for the bump 
suffered maximum number of scatterings around 8 r^. From the temperature contours, we find that the 
cloud temperature around 8 is ~ 100 keV. This explains the existence of the bump. In Fig. [T5b . we 
consider all the outgoing photons independent of their energies. The difference between the two cases is not 
so visible. This shows that the bulk velocity contributes significantly to produce the highest energy photons. 

In Fig. [T31 we explicitly showed the effects of the bulk velocity on the spectrum. We note that the bump 
disappears when the bulk velocity of the electron cloud is chosen to be zero (Curve marked 2). This fact 
shows that the region around 8 Vg in presence of the bulk motion behaves more like a black body emitter, 
which creates the bump in the spectrum. Since the photons are suffering large number of scatterings near 
this region (8 Vg), most of them emerge from the cloud with the characteristic temperature of the region. The 
effect of bulk velocity in this region is to force the photons to suffer larger number of scatterings. This bump 
vanishes for lower density cloud (low rhh) as the photons suffer lesser number of scatterings. The photons 
which are scattered close to the black hole horizon and escape without any further scattering, produce the 
high energy tail in the output spectrum. Curve 3 of Fig. [Ml shows the intensity spectrum of Case Ih (Table 
1), when there are zero bulk velocity inside 3 Vg. We find that in the absence of bulk velocity inside 3 r^, 
the high energy tail in the Curve Ih vanishes. This is thus a clear signature of the presence of bulk motion 
Comptonization near the black hole horizon. 
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Figure 12: Variation of the spectrum with the increase of the halo accretion rate, keeping the disk rate 
(rhd = 10) and angular momentum of the flow (A = 1) fixed. The case number for which a spectrum is 
drawn is marked on it. The unmarked plot is the injected spectrum. The resulting spectrum becomes harder 
for the higher values of rh/j. 





Figure 13: Number of scatterings inside the spherical shell between R to R + SR {5R ~ 1.4). The light and 
dark shaded histograms are for the cloud with and without bulk velocity, respectively, (a) Only the photons 
emerging from the cloud with energies E, where 50 keV < E < 150 keV, are considered here, (b) All the 
photons emerging from the cloud are considered here. Parameters used: rhd = 1, rfih = 10 and A = 0. 
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Figure 14: The spectrum for the Case Ih. The curves marked 2 and 3 give the spectra when the bulk velocity 
of the electron is absent for the whole cloud and for the cloud inside 3 Vg, respectively. The curve marked 1 
gives the injected spectrum. The bulk motion Comptonization of the photons inside the 3 Vg radius creates 
the hard tail. The bump near 100 keV is a combined effect of the temperature and bulk velocity of the rest 
of the cloud. 

5 Summary and Discussions 

In this paper, wc have extended our previous work using Monte Carlo simulations (Ghosh et al. 2009, 2010) 
to include the effects of Comptonization on the dynamics of the accreting halo having zero and very low 
angular momentum. We studied the properties of the emerging spectrum from the Chakrabarti-Titarchuk 
model of a two component flow, one component being the Keplerian disk on the equatorial plane and the 
other component is the low angular momentum accreting halo, which is acting as the Compton cloud. In 
Table 1, we have given the parameters of all the cases which were run. We note that as we enhance rhd, Ninj 
see, §2.2) is also enhanced, increasing the number of photons Nsc undergoing Compton scattering. 
If we keep the halo rate rhh fixed, then increasing rhd increases N^nsc, the number of photons escaping 
from the disk while keeping p almost unchanged. For Cases l(a-d) and 2(a-d), the percentage of photons 
undergoing scattering p is ~ 20% and ~ 19%, respectively. When we increase rh/,, keeping constant, this 
percentage increases rapidly due to the decrease in Nunsc- We see from Table 1, as rhh is increased from 
0.5 to 10, keeping = 1 [Cases le, la, l(f-h)] percentage of scattered photon increases from ~ 11% to 
76%. The same situation prevails for the A = 1 cases [Cases 2e, 2d, 2(f-h)], where p increases from ^ 11% 
to ^ 51%, for the increase in -rhh from 0.5 to 5, keeping -rhd = 10. N^h remains almost constant if we keep 
the halo rate constant. If we increase the halo rate, N^h increases rapidly, because increase in mu increases 
the density of the cloud and thus pushes the photons towards the black hole. As we increase rrid, the cooling 
time to decreases, since with the increase of rhd number of soft photon increases. Thus, the cloud cools down 
at a faster rate. We also notice that to decreases as rhu increases, due to the increase of Ng^- Spectral index 
a increases with the disk rate for a fixed halo rate, and it decreases with halo rate for a fixed rhd. This 
can be explained by the fact that as we increase rhd, the electron cloud becomes cooler, the spectrum gets 
softened. On the other hand, when we increase the number of hot electron inside the cloud (i.e. m/j), for 
the same Ninj we get a hotter system. This makes the spectrum harder. These results are consistent with 
the Chakrabarti-Titarchuk scenario of two component accretion. 

Our major conclusions are the followings: 
i) In the presence of an axisymmetric disk which supplies soft photons to the Compton cloud, even an 
originally spherically symmetric accreting Compton cloud becomes axisymmetric. This is because, due to 
the higher optical depth, there is a significant cooling near the axis of the intervening accreting halo between 
the disk and the axis. 
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ii) Due to the cooling effects close to the axis, the pressure drops significantly, which may change the flow 
velocity up to 25%. This effect becomes more for low angular momentum flows which produces shock waves 
close to the axis. The post-shock region cools down and the outflow falls back to the disk. This shows that 
the Chakrabarti-Manickam (2000) mechanism of the eflPects of Comptonization on outflows does take place. 

iii) The emitted spectrum is direction dependent. The spectrum along the axis shows a large soft bump, 
while the spectrum along the equatorial plane is harder. 

iv) Photons which spend more time (up to 100 ms in the case considered) inside the Compton cloud produce 
harder spectrum as they scatter more. However, if they spend too much (above 100 ms) time, they transfer 
their energies back to the cooler electrons while escaping. These results would be valuable for interpreting 
the timing properties of the radiation from black hole candidates. 
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